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We investigate transport properties of the junctions in which the graphene nanoribbon with 
the zigzag shaped edges consisting of the N legs is sandwiched by the two normal metals by 
means of recursive Green's function method. The conductance and the transmission probabil- 
ities are found to have the remarkable properties depending on the parity of N. The singular 
behaviors close to E — with E being the Fermi energy are demonstrated. The channel filtering 
is shown to occur in the case with N = even. 
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Recently, graphene-based materials with nano-meter 
sizes have been attracting much attention in both funda- 
mental and applied sciences. Among these materials, the 
graphene nanoribbon with zigzag shaped edges, which 
is abbreviated to zigzag GNR in the following, has fasci- 
nating peculiar properties as follows. 1-3 The zigzag GNR 
has a metallic band structure irrespective of the width 
N in the sense that the energy gap does not appear at 
the Fermi energy in the absence of doping, E = 0. How- 
ever, the asymptotic form of the energy dispersion near 
E = is written as E oc ±\k — ir/a,\ N with a being the 
lattice spacing and then Drude weight in the absence of 
doping becomes zero due to the vanishing Fermi veloc- 
ity for N > 2. Therefore, according to Kohn's criterion, 4 
the system without doping is classified into the insula- 
tor from the point of view for transport properties. Such 
characteristic properties are due to the fact that the one- 
particle states near E — are localized around zigzag 
edges. 

Transport properties of the zigzag GNR applied to the 
external potential have been studied theoretically, 5-8 and 
it has been found that the parity of the width N remark- 
ably affects the transport properties. In the case of TV = 
even, the zigzag GNR has the reflection symmetry in 
the transverse direction, and as a result, the parity of 
the wave function in the direction must be even or odd. 
Since only the scattering processes conserving the parity 
of the wave function are allowed under the external po- 
tential not varying in the transverse direction, the filter- 
ing of the scattering processes happens depending on the 
parity of the wave function of the transverse direction. 8 
Note that such an interesting phenomenon depending on 
the parity of the width N have been found in the persis- 
tent current of the isolated ring pierced by the magnetic 
flux. 9 However, the transport properties of the junctions 
which consist of the zigzag GNR and the normal metals 
have not been investigated. 

In the present paper, we investigate the transport 
properties of the junctions shown Fig. 1 where the zigzag 
GNR with the width N and the length 2N L is sand- 
wiched by the two ideal leads expressed by the regular 
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square lattices. The Hamiltonian is written as follows, 




Fig. 1. The present model where the zigzag GNR written by the 
thick lines is sandwiched by the ideal leads expressed by the 
regular square lattices, which is expressed by the thin lines. Here, 
the width and length are given by N and 2Nj J , respectively, and 
a is the lattice spacing of the zigzag GNR. The filled (open) 
circles in the zigzag GNR region express the A (B) sublattice. 
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where \i) is the localized state at the site i and = —t 
if i and j are nearest neighbors, otherwise ty = 0. The 
hopping in the GNR region is assumed to have the same 
value t as that in the normal metals for simplicity. The 
conductance in unit of 2e 2 /h ( — e < : electronic charge, 
h : Planck constant), g, written by Landauer formula 



(2) 



is calculated by the recursive Green's function method. 10 
Here, T^^ is the transmission probability from the inci- 
dent channel v to transmitted one fi, both of which are 
defined in the ideal leads. Note that we have confirmed 
that the conservation law J^fj, {Tn,v + R^.v} = 1 holds in 
all the numerical calculation in order to check validity of 
our results where R^, v is the reflection probability from 
the channel v to fj,. 

At first, the cases without doping, i.e., E/t = are 
discussed. 11 We show the conductance g as a function of 
the length Nl for N — even in Fig. 2 (a) and N = odd in 
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Fig. 2 (b), respectively. We can observe obvious qualita- 
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Fig. 2. The conductance g at E/t = as a function of Nl for 
N = 2, 4, 6, 8, 10 (a) and for N = 3, 5, 7, 9, 11 (b). In (a), the 
dotted lines express the asymptotic behaviors for Nl ^> 1, all of 
which are proportional to N// 2 . 
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Fig. 3. (Color online) The conductance g at E/t > as a function 
of 2k F L/ir for N = even (a) and for N = odd (b) with L = (N L + 
l/2)a being the length of the zigzag GNR region and Nl = 200. 
In each figure, the right end of horizontal axis 2kpL/n = 401 
corresponds to the E/t = case. 



tive discrepancy between N = odd and N — even in the 
asymptotic behavior for Nl 3> 1. In the case with N — 
odd, the conductances approach to unity, which indicates 
that only the one channel survives in the limit. On the 
other hand, those for N — even decrease in proportion to 
Nl 2 . Actually, we can successfully derive the expression 
g ~ 6/-/V| by obtaining the recursive Green's function 
analytically in the case of N = 2. Note that the trans- 
mission probabilities are obtained as Tx,i — T2.2 — 3/A^£ 
and Ti : 2 = Ti,\ — in this limit. 

Next, the cases with E/t 7^ are investigated. Due to 
the particle-hole symmetry, the transport properties for 
E/t > are the same as those for E/t < 0. Therefore we 
concentrate on only the E/t > cases unless explicitly 
noted. In addition, we discuss the energy regions where 
the one band gets across the Fermi energy. Here, there 
exits the one-to-one correspondence between the energy 
and the Fermi wavenumber &f ; kp is equal to 7r/a at 
E/t = and decreases with increasing the energy. There- 
fore, the conductances are investigated as a function of 
the Fermi wavenumber instead of the energy. In Fig. 3, 
we show the conductances as a function of 2k-pL/n for 
N = even (a) and N = odd (b) with the fixed length 
Nl = 200. In this analysis, we use L = (Nl + l/2)a as 
the length of the zigzag GNR region L. In each figure, the 
right end of the horizontal axis is given by 2k-pL/ir = 401 
and it corresponds to E/t = 0, i.e., kp = n/a. We can 



see the oscillating behavior in the both cases except near 
E/t = 0. The interval between one peak and the near- 
est neighbor one is given by A(kpL) = ir. Therefore, the 
oscillation is considered to be originated from the inter- 
ference between the electron waves at the Fermi energy. 
It should be noted that such oscillating behavior disap- 
pears near E/t = 0. In the cases of N=2 and 3 (4 and 
5), the peak which should appear at 2kpL/ir = 401(400) 
does not exist. With increasing the width of the zigzag 
GNR, the number of missing peaks increase. For exam- 
ple, the two peaks expected at 2kpL/ir — 399 and 401 
respectively are missing in the cases of = 6 and 7. 

Here, the transmission probabilities T M jV are discussed 
in detail for the N = even cases. In these cases, as was 
already noted, the band selective filtering is known to oc- 
cur in the transport properties of zigzag GNR under the 
external potential uniform in the transverse direction. 8 
The origin of the filtering is the reflection symmetry in 
the transverse direction, which makes the wave function 
in the direction either even or odd function. Since the 
present external leads has also reflection symmetry, the 
wave function of the channel defined in the external leads 
is either the even or odd function. Therefore, the filter- 
ing depending on the channel indices is expected for the 
N = even cases. We show the transmission probability 
T^ v for N = 4 as a function of Nl for E/t = 0.1 (a), 
E/t = -0.1 (b) and E/t = 0.0 (c) in Fig. 4. In the case 
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Fig. 4. (Color online) The transmission probability Tn, >v together 
with the conductance g for N = 4 as a function of in the 
cases of E/t = 0.1 (a), E/t = -0.1 (b) and E/t = 0.0 (c). 

of E/t = 0.1, the transmission probabilities T oc id,odd be- 
comes finite for Nl 3> 1, whereas only T evonjeven survives 
for E/t = —0.1. The other transmission processes such 
as those from the even channel to the odd one and vice 
versa are suppressed in the both cases. These results can 
be understood in terms of the parity of the wave func- 
tion in the transverse direction; The wave function of the 
channel //, it M (7) (I = 1, 2, • • • , N), in the external leads 
is given as follows, 

= \j N \ 1 sinky(fi)l, (3) 

and it's parity is positive for fj, = odd, whereas negative 
for fi = even. On the other hand, in the zigzag GNR 
with the width N = even, among the two bands located 
around E/t = the wave function for E/t > is the 
even function in the transverse direction, whereas that 



Fig. 5. (Color online) The transmission probability T^^ together 
with the conductance g for N = 4 and Nl = 200 as a function 
of 2k-pL/n. The right end of the horizontal axis 2k^L/n = 401 
corresponds to E/t = 0. 

for E/t < is the odd one. 8 Since the parity of the wave 
function in the transverse direction should be conserved 
when passing through the zigzag GNR, the channels with 
positive (negative) parity can pass through for E/t > 
(E/t < 0). At E/t = 0.0, the two bands touch with each 
other and two states are degenerate. Therefore, the two 
kinds of channels, i.e., from even to even and from odd 
to odd, are alive in the case of E/t = 0.0. Actually, the 
analytical calculation for N = 2, T ls = T 2i2 ^ 3/iV| 
supports the conclusion shown above. 

Now, we investigate this channel filtering in detail close 
to E/t = where, as has been already discussed above, 
the anomalous phenomenon that the oscillation of g as 
a function of kpL disappears is observed. In Fig. 5, the 
transmission probabilities for N — 4 and Nl — 200 are 
shown as a function of 2kpL/n as well as the conduc- 
tance g. Here, the E/t ^ case is discussed and the right 
end of the horizontal axis corresponds to E/t = 0.0. The 
last peak in the conductance oscillation is observed at 
2fcpi/7r = 398 and the peak which should be observed 
at 2k^L/-K = 400 diminishes. In the region, we can see 
that the transmission processes from the even channel to 
the even one, which is prohibited due to parity conserva- 
tion, become possible. The regions where the oscillation 
disappears correspond to those with the breakdown of 
the parity conservation. Therefore, the two phenomena 
can be considered to come from the same origin. 

Thouless discussed that the energy levels of the finite 
system attached to the leads are shifted by the amount 
of H/t, where r is the time it takes for an electron to 
move to the end of the system. 12 The conductance can 
be understand from the point of view that the states 
within the shift can pass through the sample. 13 Though 
he discussed the case with random impurities where the 
motion of the electron is diffusive, such a shift of the en- 
ergy levels in the GNR region of the present system can 
be expected to occur in the presence of the ideal leads. 
As a result, the anomalous phenomena seen in E/t = 
can be expanded to the E/t ^ region. However, at 
present, it is unclear why the conductance at E/t ~ 
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shows the discrepancy in the Nl dependences between 
N = even and N — odd case. It is likely that the dis- 
crepancy in the transport properties originates from that 
in coupling between the GNR and the normal metal be- 
cause the qualitative difference between N = even and 
odd cannot be observed in the band structure of the iso- 
lated GNR. Further investigation on the electronic states 
of GNR attached with the ideal leads is necessary in or- 
der to understand the transport properties of the normal 
metal - zigzag GNR - normal metal junctions. 

In conclusion, wc investigated the transport proper- 
ties of the junction where the zigzag GNR with width 
N and length 2Nl is sandwiched by the normal metals 
expressed by the regular square lattices. The transport 
property at E/t — shows different asymptotic behavior 
as the a function of Nl, i-e., for Nl 3> 1 the conductance 
with TV = odd approaches to unity whereas that in the 
case of N = even decreases in proportion to N£ 2 . The 
conductance at E/t ^ shows oscillating behavior as 
a function of the energy, which is due to the interfer- 
ence of the electron wave at the Fermi energy. Such a 
oscillation disappears close to E/t = 0. In addition, the 
channel filtering occur in the case of N = even with the 
reflection symmetry in the transverse direction. The fil- 
tering is found also to show unexpected behavior close 
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to E/t = 0. 
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